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Abstract 



We extend and test empirically the multifractal model of asset returns based on a multi- 
plicative cascade of volatilities from large to small time scales. Inspired by an analogy between 
price dynamics and hydrodynamic turbulence [Ghashghaie et al, 1996; Arneodo et al, 1998a], 
it models the time scale dependence of the probability distribution of returns in terms of a 
superposition of Gaussian laws, with a log-normal distribution of the Gaussian variances. This 
multifractal description of asset fluctuations is generalized into a multivariate framework to ac- 
count simultaneously for correlations across times scales and between a basket of assets. The 
reported empirical results show that this extension is pertinent for financial modelling. Two 
sources of departure from normality are discussed: at large time scales, the distinction between 
discretely and continuously discounted returns lead to the usual log-normal deviation from nor- 
mality; at small time scales, the multiplicative cascade process leads to multifractality and 
strong deviations from normality. By perturbation expansions, we are able to quantify precisely 
on the cumulants of the distribution of returns, the interplay and crossover between these two 
mechanisms. The second part of the paper applies this theory to portfolio optimisation. Our 
multi-scale description allows us to characterize the portfolio return distribution at all time 
scales simultaneously. The portfolio composition is predicted to change with the investment 
time horizon (i.e., the time scale) in a way that can be fully determined once an adequate mea- 
sure of risk is chosen. We discuss the use of the fourth-order cumulant and of utility functions. 
While the portfolio volatility can be optimized in some cases for all time horizons, the kurtosis 
and higher normalized cumulants cannot be simultaneously optimized. For a fixed investment 
horizon, we study in details the influence of the number of periods, i.e., of the number of rebal- 
ancing of the portfolio. For the large risks quantified by the cumulants of order larger than two, 
the number of periods has a non-trivial influence, in contrast with Tobin's result valid in the 
mean-variance framework. This theory provides a fundamental framework for the conflicting 
optimization involved in the different time horizons and quantifles systematically the trade-offs 
for an optimal inter-temporal portfolio optimization. 



*We are grateful to E. Bacry, U. Frisch and L. Martellini for helpful discussions. 



1 Introduction 



Inspired by an analogy with turbulent cascades in hydrodynamics comparing the energy flow to an 
information transfer, recent empirical works have shown that return volatilities exhibit long-range 
correlations organized in a hierarchical way, from large time scales to small time scales [Ghashghaie 
et al, 1996; Arneodo et al, 1998a; Muzy et al, 2000; Breymann et al., 2000]. The Olsen group 
in Zurich discovered independently this phenomenon which they called the HARCH effect [Miiller 
et al., 1997; Dacorogna et al, 1998]: in the foreign exchange market, the coarse-grained volatility 
predicts the fine-grained volatility better than the other way around. They also found this effect 
for the implied forward rates derived from Eurofutures contracts [Ballocchi et al, 1999]. 

The underlying cascade or hierarchical structure provides also a natural explanation and a model 

for the multifractal description of stock market prices documented by several authors [Fisher et al, 
1997; Mandelbrot, 1997; Vandewalle and Ausloos, 1998; Brachet et al, 1999; Bershadskii, 1999; 
Ivanova and Ausloos, 1999; Mandelbrot, 1999; Schmitt et al. 1999; Pasquini and Serva, 2000] (see 
however Bouchaud et al. (2000)). 

The first purpose of the present paper is to provide additional empirical confirmations that such 
hierarchical description accounts very well for the return statistics, especially for the strong lep- 
tokurticity of the probability density functions at small scales and for the volatility serial correla- 
tions. This hierarchical framework shares a crucial advantage with models based on the standard 
geometrical Brownian motions or on Levy-stable processes, in that it remains a "scale-free" descrip- 
tion which can thus be applied to any time scale. As we shall see, this "self-similarity" property 
is interesting for the problem of portfolio optimization, in particular for the multi-period selection 
problem. 

The portfolio optimization problem consists in finding the optimal diversification on a set of possibly 
dependent assets in order to maximize return and minimize risk. In its simplest version, one assumes 
a single period horizon for all investors which, together with the hypothesis that returns are normally 
distributed, leads to the standard Markowitz's solution [Levy and Markowitz, 1979; Kroll et al, 
1984]. The first results in the theory of optimal multi-period portfolio selection showed that, within 
the Gaussian hypothesis of return distributions, the investor optimal sequence of portfolio through 
time is stationary, with constant proportionate holdings of each included asset leading to a constant 
expected return and risk per unit invested wealth. In other words, the optimal asset weights in the 
portfolio are time-scale independent. This theorem has been criticised as being not true in general 
[Stevens, 1972]. 

Many studies have investigated the impact of time horizon on the portfolio selection. Allowing 
for investors with different planning horizons, Gressis et al. [1976] have shown that, in absence of 
the riskless security into the portfolio, the single-period efficient frontier also provides K periods 
efficient frontiers. However, Gressis et al. [1976] find that, with the inclusion of the riskless security, 
the equivalence among efficient frontiers for different horizons is no longer valid. Departing from 
the Gaussian model, Arditti and Levy [1975] generalized the mean-variance approach to a three- 
moment efficiency analysis and showed that all one-period portfolios are not necessarily multi-period 
efhcient because, even if a stock has a symmetric single-period return distribution, its multi-period 
distribution may be highly skewed. A considerable body of research has explored how portfolio 
composition depends on the investment horizon. It has been shown that portfolio composition 
indeed depends on the investment horizon and that any simple characterization of the relationship 
is treacherous [Gunthorpe and Levy, 1994; Ferguson and Simaan, 1996]. Using a choice criterion 
that is consistent with the traditional utility approach but which is more amenable to a multiperiod 
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environment, Marshall [1994] has shown that investors should choose progressively less risky single- 
period portfolios as their investment horizons shorten, even if they do not become more risk averse, 
both in the presence and in the absence of a riskless asset. Tang [1995] examined the effect of 
investment horizon on international portfolio diversification using stock indexes of 11 countries. 
He finds that correlation coefficients between stock indexes increase in general with an increase in 
the investment horizon, suggesting that diversification benefits are reduced. The fact that various 
stock markets are more correlated over a longer investment horizon may imply that different stock 
markets adjust to each other with a delayed pattern. Bierman [1997; 1998] has investigated whether 
the risk of a stock portfolio increases or decreases as the investment horizon lengthens and whether 
the portfolio mix depends on the horizon. By moving beyond the logarithmic utility function and 
by recognizing that different investors have different risk preferences, Bierman [1998] finds that it 
is possible for stocks to be risky in each time period and to reduce risk by increasing the length 
of the investment horizon. He also finds that stocks can be rejected if they are to be held for one 
time period, but can then be accepted if they are to be held for more than one time period. At 
the origin of these results is the fact that "discrete" returns r calculated as the ratio of the price 
difference over a given time interval r over the price are only approximately proportional to r times 
the continuously discounted return rj. The difference is negligible at short times r, corresponding 
to the validity of the expansion of e^'^ = 1 + rjT. At large time scales such that ijt is no more 
small compared to 1, the second-order term in the expansion becomes important and effets of 
non-normality in the distribution of r become important. In the present paper, we study another 
additional cause for non-normality, whose strength grows as short time scales rather than at large 
time scales. 

The paper is organized as follows. In section 2, we recall and make precise the statistical description 
of the fluctuations of the returns of a single asset in terms of the multiplicative cascade model. The 
few parameters involved in this model are then interpreted and estimated for a set of high frequency 
time series. Then, we propose an extension of this framework to vector valued, i.e. multivariate 
processes corresponding to a basket of assets in a portfolio. This can be done in a very natural 
way and simple statistical tests are proposed to check for its relevance. In section 3, we address 
the portfolio problem for such a set of multifractal assets. We first discuss the simplest case of 
uncorrelated assets where the statistical properties of the portfolio returns can be estimated by 
means of a cumulant expansion. The problem of portfolio optimization is studied in section 4 using 
the utility function approach and the higher-order cumulants of the distribution of portfolio returns. 
A generalisation of the efficient frontier is proposed and studied as a function of time-scales. Section 
5 presents the application of the correlated multivariate multifractal model to the portfolio theory. 
Conclusion and prospects are given in section 6. 

2 Multifractal description of asset returns 

In the following, we show how multifractal statistics can be described in a simple way using only a 
few parameters. This description has been developed both for finance data [Arneodo et al, 1998a] 
and in the field of fully-developed hydrodynamic turbulence [Arneodo et al, 1998b; 1998c, 1999]. 
As mentionned above, the similarity between turbulence and finance has been suggested by some 
recent studies [Ghashghaie et al, 1996; Arneodo et al, 1998a; Muzy et al, 2000] and relies on the 
concept of multiplicative random cascades: the fluctuation at some fine time scale ti is the product 
of the fiuctuation at a coarser scale T2 by a random factor whose law depends only on the scale 
ratio T2/T1. This picture leads to the integral description proposed by B. Castaing and collaborators 
[Castaing et al., 1990; Chabaud et al, 1994] which links the probability density function (pdf) of 
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fluctuations at scale ri to the pdf of fluctuations at scale T2. Let us now introduce this formalism 
for the pdf's of asset returns for all time scales. 



2.1 Formulation of the multifractal model of the returns of a single asset 



Let Pi{t) be the price of asset i at time t, where time is counted for trading days in multiples of a 
fundamental unit (say days). With the notation ryj(f,r) = In (^^j^^ttj); the return rj(t,r) between 
time t — T and t of asset i is defined as: 

r.(t,r) = ^i(i-^i(if^ = e-(*-)-l, (1) 

Pi\t T) 

We thus distinguish between the "continuous return" r]i{t, r) and the "discrete" return ri{t, r). The 
difference between them is essentially captured by the second order correction |[r/j(t, r)]^, which 
becomes non-negligible only for large time scales (see below). 

We describe the distribution of the variable r]i{t,T) by using the representation of Castaing et 
al. [1990; Chabaud et al, 1994] inspired by the above mentionned cascade picture for hydrody- 
namic turbulence. This distribution is represented, for time scale ri, as a weighted sum of dilated 
distributions of the r/'s at the intermediate coarser scale T2 (ri < T2): 

/+00 
du e"" Gi^r2^ri{u) Pi,T2(e~"(??i - fJ.iT2)) , (2) 
-00 

where //jr is the mean value of rji at scale r (it is supposed to be linear with r consistently with 
experimental observations) and Gi^r2^Ti{u) is the self- similarity kernel whose shape is assumed to 
depend only on the ratio T2/T1. As we shall see below, the exponentials e~" are just convenient ways 
of simplifying the parameterization of the kernel. This model (|2|) can be derived from a cascade 
model according to which ri(t,Ti) is written as 

rjit,n)=ait,n)X{t) , (3) 

where X{t) is a ri-independent Gaussian random variable and the "stochastic volatility" a{t,Ti) is 
described by a multiplicative cascade [Arneodo et al, 1998a]: 

a{t,Tl) = Wr2^r^(Tt,T2 ■ (4) 

It is a well-established observation that, at sufficiently large scale, distributions of r/ become Gaus- 
sian [Campbell et al., 1997]. Let Tj be such a large scale for which Pi^TiiVi ~ l^iPi) is a Gaussian 
AA(0,cr|^). Then, for t <Ti, expression (|^) can be rewritten as 

/ + 00 
du Gi^T,^-r{u) Mi^T^iVi - l^iTi-.e^'^a^^) . (5) 

Pi,T{'n " l^'i'^) is thus a weighted sum of Gaussian distributions and the kernel Gi^Ti-^riu) can be 
identified with the pdf of the logarithm of the multiplicative weights u = IniyVTi^r) introduced in 
Eq.(|^, the standard deviation at scale r being: 

(Ji = cVt, . (6) 

Let /ii(r) and A?(r) be respectively the mean value and the variance of Gj.-r-^r- We consider the 
simplest possible case, specifically where Gi^Ti^rW) is itself Gaussian: 
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From the semi-group property resulting from the cascade equation, WVa-^n = Wr3^T2Wr2^Ti , and 
using the fact that Gi^r2^Ti depends only on the ratio T2/T1, it is easy to show that /ij(r) and A?(r) 
are necessarily linear functions of ln(T): 

/ii(r) = ai + hiluT , (8) 
[A,(r)]2 = b.-X^lnr, (9) 

where ai,bi,Xi and hi are independent of r. Tj is the coarse scale defined by a vanishing variance 
XfiTi) = 0. This time scale Tj is the "integral" scale at which the cascading process of the volatility 
begins. Performing the change of variable u — > u — hi{Ti) amounts to keep Eq.(|5|) unchanged and 
to redefine Eqs (|8|) and (P) as : 

hi{T) = hiH-r/Ti) , (10) 
[X,{T)f = -\j\n{r/T,). (11) 

As shown in the next section, these equations account very well for empirical observations. 

Let us note that this description (^) with (^) and ( pT| ) contains the standard Gaussian random walk 
model of stock market returns as the special case : 

hi = \. (12) 
Ai = 0. (13) 

Indeed, for Aj — > 0, Gi^T^^riu) becomes a delta-function centered on u = hi, i.e. the volatility @ 

at scale r is o"i = o"T,e" = o"T,e ' = gt^ ( ^ 1 , which recovers the standard square-root diffusion 
law. In this case, the pdf of asset "continuous" returns Pi^ri'ni ~ f^iT) given by Eq.(|5|) is Gaussian. 

2.2 Empirical tests 

In order to test the previous formalism empirically, we need to define estimators for the parameters 
Tj, 0"!^ , hi and Xf. For this purpose, the so-called "structure functions" or g-order moments {r]i{t, r)) 
of the "continuous" returns are very useful quantities (the brackets define the average with respect 
to a running window along the time series). For the purpose of analyzing data with multiple scales, 
the dependence of {r]1(t,T)) as a function of time scale r for various q's is very rich in information 
[Frisch, 1995]. For the hierarchical model (^) with the Gaussian ansatz given by Eq.(^), the gth- 
order structure function of the centered variable r]^{t, r) = r]i{t, r) — /i^r (whose law is denoted using 
prime) 

f+00 

M,{q,T) = {rj't)= / rl^PlM)d^,. (14) 



can be calculated analytically with the change of variable R = rj^e 

/ i'+oo \ / r+co \ 

Miiq,T)=i^J_ e'^^G,,T,^riu)duj i^J^ R'^PlT,iR)dRj . (15) 



The second factor in the r.h.s. of Eq.(|l5|) is the gth-moment of the distribution at the integral 
scale. All the dependance on the time scale r is found in the first factor that can be computed 
using Eqs and (0): 

+00 / T- \ Q{g) 



00 



e'^'^G,,T,^riu)du =(-] , (16) 
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with 



C^{q) = Kq-^q' • (17) 



We thus get 

M,{q,T)='yj_^WPlrr^{R)dRj (^-j . (18) 

The dependence of the gth-order moments of "continuous" returns as a function of time scale is 
a pure power law, resulting from the logarithmic behavior as a function of time scale of both the 
mean and the variance of the Gaussian self-similarity kernel (Eqs. ( [To| ) and ([Tl|)). The exponent 
C,i{q) is related to the cumulant generating function of the kernel, as shown by Eq.(^). In general, 
the centered moments of odd-orders at the integral scale are vanishing due to the approximate 
symmetric structure of the pdf's. In particular, the first-order moment, the centered "continuous" 
return, is zero by definition. From expression (|T8|), the same property follows for {rj'l). In order to 
measure the power law dependence given by the second term (r/Ti)'^^'^^ of the r.h.s. of Eq.(18), it 
is convenient to take the absolute value of the "continuous" return and calculate : 

• (19) 



In order to determine the two parameters hi and Aj, we measure the "multifractal spectrum" Ciiq) 
by performing a linear fit of In Mj (q, r) as a function of In r, for a wide range of values q. The slope 
of the linear fit gives the exponent Ciiq)- The slope at the origin of Ciiq) as a function of q provides 
an estimation of the parameter hi. The fit of Qiq) ~ hiq then provides a check for the importance 
of the quadratic correction proportional to A^. Recall that, for a Gaussian pdf of "continuous" 
returns, A? = and the exponents Ci{q) linear in q. The quadratic dependence, often called 
"multifractal" in the literature [Prisch, 1995], is a signature of the multi-scale structure. 

In Fig. m, we report such an analysis for both the S&PSOO index future and Japanese yen/US dollar 
exchange rate time series shown at the top in the two panels (a). The original intraday series have 
been sampled at a 10 mn rate and seasonal effects on the volatility have been removed. One can 
see in Figs |l](b) that the moments Mj(g, r) do exhibit the predicted scaling (Eq. (p!9[)) over a large 
range of scales, of which three decades are here shown. The multifractal nature of these series is 
illustrated in Figs ||(c) which shows the ratio ^^^i^-^q for g = 2, 3, 4 and 5 in log-log coordinates. 

For linear C,i{q) functions, the ratios j!^^^i^^q should be constants as a function of time scales. A 
departure from a constant thus qualifies a multifractal behavior. The dependence of the exponent 
Clio) reported in Figs ||(d) is very well fitted by the parabolic shape given by Eq. (|T7|). We find 
/ij — A? = 1/2 to a very good approximation for all assets which, according to Eq. (p!7|), ensures 
that the volatility or variance {r]'1) ~ t as for the standard geometrical Brownian model, i.e. does 
not exhibit an anomalous scaling. 

In order to test for the existence of a cascading process and to estimate the integral time Tj, we 
refer to [Arneodo et al, 1998a] where it is shown that, for a cascade process, the covariance of the 
logarithms of the r^j's at all scales should decrease as a logarithmic function. More precisely, if one 
defines the magnitude of asset i at scale r and time t as LOi{t,T) = ln(|r/j(t, r)|), then one should 
have, for At > r: 

Ci{At, r) = Coviujiit + At, r), c^i(t, r)) ~ -Af In(^) . (20) 

-Li 

This behavior is checked for the S&P500 and JPY/USD series in Fig. |2| where we have plotted 
C(At,r) versus In(At) with r = lOmn. The linearity of these plots is reasonably well verified. 
The values of A? and In(Tj) can be obtained respectively from the slope and the intercept of these 
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Series 




T (year) 


/"i (year ^) 


o"r, (year ^) 


German Government Bond (F) 


0.028 


1.2 


-1.6 10"^ 


1.7 10"^ 


FT-SE 100 Index (F) 


0.017 


1.1 


-3.3 10"^ 


7.6 10-5 


Japanese Yen (F) 


0.032 


1.0 


-1.1 10-^ 


3.8 10-^ 


S&P 500 (F) 


0.018 


4.5 


5.7 10--^ 


1.8 10-4 


Japanese Government Bond (F) 


0.121 


1.3 


9.2 10-^ 


4.0 10-^ 


Nikkei 225 (F) 


0.030 


1.5 


5.5 10-^ 


3.9 10-^ 


French CAC40 (mean) (S) 


0.020 


1.9 


1.2 10-1 


5.6 10-^ 


Dow-Jones (mean) (S) 


0.010 


1.7 


1.7 10-1 


4.6 10"^ 



Table 1: Estimates of the multifractal parameters for intraday (S&P 500 and Japanese Yen futures) 
and daily time series. The values measured at different time scales are consistent with each other. 

straight lines. In Table 1, we have reported the estimates of the parameters A?, al^, , Hi and Tj for 
a set of high frequency future time series and the mean values of those parameters for daily time 
series of stocks composing the Dow Jones Industrial Average index and the french CAC40 index. 
Let us note that the errors on the estimates of Tj can be very large since we get only one estimator 
of In(Tj). For all series, we have checked that the relationship hi — Xf = 1/2 is reasonable so we 
have not reported the estimated values of hi. We can remark that the values of A? are relatively 
close to 0.02 and the values of Tj are in the range 1 — 2 years. These values can thus be considered 
as representative of market multifractality. 

2.3 The multivariate multifractal model 

Up to now, our emphasis has been on the correlation structure within each asset separately. We have 
not investigated the impact of correlations across assets and have quantified only the "diagonal" 
multifractal features as it is the case for uncorrelated assets. Since one of the natural applications 
of the characterization of distributions of returns is the quantification and selection of portfolios 
made of a basket of assets, it is important to generalize our framework to account for the possible 
existence of inter-asset correlations, in addition to the multi-scale time correlations. As it is well- 
known, optimizing a portfolio relies on two effects, the law of large numbers and the possibility 
to counter-act the negative effect of one asset by using (anti-) correlations. Towards a practical 
implementation of our model, it is thus essential to offer a generalization of the multifractal cascade 
model that accounts for the correlations across assets. 

When dealing with Gaussian processes, the lack of independence is entirely embodied in the corre- 
lation function. Along this line of thought, we now propose a phenomenological model that extends 
the previous one to the case of correlated assets. We would like to generalize Eqs. (|5|), (0), (|8|) and 
(^ to multivariate distributions. In this purpose, let us write *f/ the line vector (r/i, ...,r/jv), where 
the subscript * stands for 'transpose', and we suppose that there exists a time scale T for which 
the multivariate distribution MN,Ti'n, Sq) is Gaussian with covariance matrix Sq. The multivariate 
Gaussian mixture that generalizes Eq.(|5|) can be written as 

Pr{f]-Tfl) = J d^'u Gr{u) Mn,t (r/ - r/2,D(5)5:oD(n)) , (21) 

where D(u) is the diagonal dilation matrix Djj = e^^6ij and Gr{u) is some multivariate probability 
distribution function. This description corresponds to a multivariate cascading process where the 
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multiplicative weights W^^tj associated with different assets can be correlated. This form (^TJ) is 
the multivariate generalization of (|5|). 

As in the previous section, let us make a Gaussian ansatz for the joint law G of the logarithm of 
these weights, i.e., 

Gr{u)=AfN{u-h{r),A{T)) , (22) 
where the mean vector /i(r) and the covariance matrix A(r) are supposed to behave as: 

hi{T) = ai + hiln{T), (23) 
Au(^) = -A?,ln(r/r,,), (24) 

with the "integral times" Tij defined by 

A,,{t = T,,) = . (25) 

Within the cascade model, the times Tij are the time scales at which some cascading process begins. 
If one considers two assets i and j, it seems natural to assume that only two situations occur: 

i/ the cascade on i and j has the same "economic" origin, and then Ta = Tjj = Tij] 
ii/ the cascades on i and j are independent and thus Ajj(r) = 0, but, a priori, Ta / Tjj. 

The returns rji can be thus organized by "classes" in such a way that the matrix A is a block 
diagonal matrix: 

L 

A = 0Az, (26) 

1=1 

where L is the number of indenpendent subspaces (cascades) of dimension Ni {J2i=i — At 
fixed I, one can always redefine "Sq such that Eqs. (p3|) and ( p^ ) become : 



h{T) = hi\n{T/Ti) (27) 



Aij(r) = -A?- \n{T/Ti) for i, j in the bloc I . (28) 



To demonstrate that this multifractal multivariate model is compatible with the empirical obser- 
vations, we first define an estimator of the elements of the covariance matrix A. Let us denote 
u}i{t,T) the magnitude, i.e., the logarithm of the "continuous" return tj for asset i at scale r: 
<^i{t,T) = ln(|ryj(t, r)|). In the multivariate framework, the structure function at scale r can be 
extended to order p, where p is a vector of N values, in the following way: 

(|,?i(t,r)r ...|w(t,r)r) = (e-(*'-)-P) , (29) 

where a;(t, r) is the vector constructed from the Wj's at scale r. Using Eq. (21), the behavior of 
the structure function as a function of r can be estimated: 



where G is the kernel quantifying the cascade from the integral scales to scale r. Using our 
assumption that G is normal with mean /i(r) and covariance A(r), the characteristic function of 
G, defined as the last term in the previous expression, is 
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From the logarithmic behavior in Eqs (^) and (p^, one deduces: 

(e'^(*'-)-P) ~ T<^P) , (30) 

with 



- 1* 1 
C{P) =h.p+- pAp = Y,Pihi - ^ ^PiPj>Hj ■ (31) 



The constant K^; depends on the values of the integral scales Tjj, /ij, Xij and on the multivariate 
distribution at the large integral scales. This power law ( |30D can be explicitely used to estimate 
the values A?^ of the elements of the covariance matrix defined by the kernel G, by using different 
set of p's. 



If one sets Pk = ^ for k ^ and pi = pj = e ^ in the expression (30), one can show that the 
covariance between the magnitudes uji{t,T) and LUj{t,T) simply behaves as: 

CoY{uJi{t,T),ujj{t,T)) ~ -A2. ln(r) + , (32) 

where Kij depends on Tij, hi and Xij. 

Alternatively, one can test the multivariate cascade ansatz from the behavior of the lagged magni- 
tude covariance: in full analogy with the monovariate case, it should behave as: 

Cij{At) = Cov{LOi{t + At, T)ujj{t, r)) = -Xfj ln{At/Tij) , (33) 

for At > T. This equation provides other simple estimators of both Xij and Tij. 

We were not able to use intraday data to test these predictions because the formalism requires 
that the different time series should be sampled at the same times. This can be alleviated in 
the future by a suitable pre-treatment that ensures the coincidence of the sampling times. Here, 
we present results for daily returns on the magnitude covariance functions for stocks taken from 
the French CAC40 index (over the period from 1992 to 1999). In Fig. |3|(a), we report in dotted 
line the typical behavior of Cij{At) versus In(At) for two assets of the CAC40 index (CCF and 
MICHELIN). Despite the relatively poor statistical convergence, one clearly sees a slow decay which 
is compatible with the logarithmic law (|33|). From the slope and the intercept of this curve, one 
gets an estimate of Xfj and Tij for these two assets. In solid line, we have plotted the average of Cij 
performed over all the pairs of assets in the CAC40 index. The logarithmic decay of the magnitude 
covariance is thus a remarkably stable feature among all pairs of assets constituting the CAC40. In 
Fig. §(b), we show the histogram of the values of Xfj estimated for 253 pairs of stocks belonging to 
the CAC40 index. Even though our estimates exhibit a large dispersion, one clearly sees that the 
distribution of Xfj is centered around A?j = 0.02. This value, which is very close to the estimates of 
A? determined in the previous section with the mono-asset analysis, can be considered as typical of 
financial stocks. This is confirmed in Fig. ^ which shows the estimates of Xfj versus the estimates 
of Tij for pairs of stocks taken both from the CAC40 and the Dow-Jones indices. Moreover, one 
sees in Fig. Q that the integral time scales Tij are clustered around T = 1 — 2 years. 



3 Characterization of the distribution of portfoho returns 

Consider a portfolio with rii shares of asset i whose initial wealth is 

N 

WiO) = Y.niPi{0) . (34) 

i=l 
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A time r later, the wealth has become W{t) = X^i^i f^iPii''') and the wealth variation is 



S^W = W{t) - W{0) = ^ nMO)^^^^^^;^j£^ = WiO) ^^ {^"^'^^^ " l) , (35) 

where 



Wi = ^ — (36) 



Ef=injPjiO) 

is the fraction in capital invested in the ith asset at time 0. Using the definition (0), this justifies 
to write the return Rr of the portfolio over a time interval r as the weighted sum of the returns 
rj(r) of the assets i = 1, ...,N over the time interval r 



3.1 Cumulants of the variables rji 

The portfolio is completely characterized by the distribution of its returns Rr for all possible time 
scales r. When the assets are assumed to be independent, expression ( pTj ) shows that the pdf of R-^ 
is obtained by the convolution of the pdf's of the returns rj(T) of the individual assets. A standard 
strategy is to calculate the characteristic functions Pi{k) of the pdf's of the returns ri{T), from 
which one extracts the cumulants and use their additive properties to get the cumulants of the 
distribution of Rr- In this section 3 and in section 4, we use the model of sections 2.1 and 2.2 of 
independent multifractal assets and turn in section 5 to the case of correlated multifractal assets 
described by the formalism of section 2.3. Let us thus estimate the cumulants of each asset return 
rj(r) at scale r. 

In that purpose, let us compute the cumulants of the "continuous" returns r/j's. The characteristic 
function of r/j is obtained from Eq.(^) with Eq.(|^ : 



oo 



From this equation, we see that the moment m\{T) of order n of the centered variable r]'^{t,T) = 
r]i{t, t) — fiiT can be easily computed 

U=o = <(r.) / du e'^" ^ , ^ exp \ ' (39) 



g^„,.=u V2^A.(r) "V 2A2(r) 

where ml^(Ti) is the n-th moment of the Gaussian law at the integral scale. The integrals in ( ^8[ ) 
and (|39| ) are proportional to the Laplace transform of the normal distribution. 

Using Eqs. (0) and (^), we obtain the following expressions for the moments of the centered 
variables rj'^: 



mUr) = (2n-l)!!4r [^J ' ^^^^ 

m*2n+i(T) = . (41) 
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Expression Eq.(pOD is the same as Eq.(^) with ( p^ ) for q = 2n. We can then easily retrieve the 
moments of the (non-centered) variables 7]i{t,T) by using the definition of r]' and inverting it as 
Vi{t, t) = r]'-{t, r) + HiT. Thus, 



for integer orders q. 



q\ 



(42) 



Note that the knowledge of the moments do not allow a unique reconstruction of the pdf. This is 
related to the fact that the characteristic function of the lognormal distribution cannot be expanded 
in a Taylor series based on the moments, because the moments grow too fast (as exp in^)) and lead 
to a diverging series. Holgate [1989] has shown how to make sense of such expansion by either 
using a finite Taylor series or resumming formally the divergent moment expansion. Using the 
entire function or moment constant method and Hardy's formula, one obtains [Holgate, 1989] 



2^— -'^T. exp 

j=0 



2J 



(43) 



where L is the Laplace transform of the lognormal distribution, which is taken at exponentially 
increasingly spaced points. Its fast decay ensures the convergence of the series against the rapidly 
diverging moments and a unique characterization of the pdf. 



The first six cumulants are linked to the moments as [Ord, 1994]: 



Cl 


= mi 


C2 


= m2 


C3 


= rns 


C4 


= 1714 


C5 


= rn5 


C6 


= ruQ 



3m2mi + 2mi , 



1 ) 



5m4mi — Wmsm2 + 20m3m^ + ciOrn^mi — 60m2m'l + 24m^ 
6m^mi — 15?n,4m2 + 30?n,4?n,^ — 10m| + \2Qm^m2mi 



-270mn 



(44) 
(45) 
(46) 
(47) 

1 , (48) 
I20m^m\ + mml 

(49) 



Using Eqs (H), (H) and (||), we then determine the six first cumulants of Notice that the 
cumulants of order larger than 1 are invariant with respect to the translation r]i{t, r) = ri[[t, r) + /i^r 
and the centered moments (pOD, ( pD can thus be used directly in expressions 



4(t 



2 



0, 

34^ 
0, 

15o-f.. 



2(h,-\1) 



6(/i,-3A2) 



4A2 



8Af 



+ 2 



12Af 



(50) 
(51) 
(52) 
(53) 
(54) 
(55) 



Expression ( |50[) retrieves the linear dependence of the mean value as a function of the time interval 
r. This result is independent of the model. Expression ( |5l| ) shows a more interesting structure: 
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for /ij — A? = 1/2 we retrieve also a linear dependence of the variance C2(t) as a function of time 
scale T. This linear dependence occurs in absence of or for weak correlations of the returns. As 
this is a ubiquitous property confirmed by all empirical studies [Campbell et ah, 1997; Mantegna 
and Stanley, 2000], this leads to the empirical constraint 

hi - A? = 1/2, (56) 

which is remarkably well verified (see Fig. ||). 

We note that, by construction of the superposition (^) of Gaussian pdf's, all odd-order cumulants 
(except the first cumulant corresponding to the average "continuous" return) are zero. The present 
hierarchical model cannot account for possible skewness in the distribution of the "continuous" 
returns rj. But as shown in the next section, this does not implies that the skewness of return 
distribution is zero. A small value of the skewness is indeed observed in empirical studies of data 
sets [Campbeh et al, 1997]. 

If one assumes that, for all assets hi ~ 0.5, then the two sources of statistical non-normality for r/, 
i.e. non zero high order cumulants, are according to our model the intermittency parameter A, and 
the integral time Tj: "quasi-Gaussian" statistics is recovered if the scale r is very close to Tj or if 
the intermittency parameter Aj is small enough so that the fluctuations of the variances at different 
scales are small. Conversely, the distribution of "continuous" returns exhibits heavy tails if Aj is 
large or if the time scale r is small compared to integral scale Tj. 

3.2 Portfolio cumulants: perturbation expansion and multifractal corrections 
to the log-normal portfolio 

In the standard portfolio theory [Markovitz, 1959; Merton, 1990], asset price time series are repre- 
sented by geometric Brownian motions, i.e the statistics of "continuous" returns is exactly Gaussian 
while that of the "discrete" returns can be approximated by Gaussian distributions by means of 
an expansion where the small parameter is time measured in year. This results from the definition 
showing that the difference between ri{t,T) and rii{t,T) is equal to |[?7j(t, r)]^ to leading order. 
Since the typical magnitude of [r]i{t, r)]^ is given by the variance of T]i(t, r) and since the mean and 
the variance of r]i{t, r) are both of the order O.lr (see stock examples in Table 1) when the time r is 
measured in year, rj(t, r) and r]i{t, r) are indistinguishable in practice for times r of the order of or 
smaller than 1 year. An alternative way of saying the same thing is that the lognormal distribution 
of rj(t, r) reduces to its approximate Gaussian representation for times less than or of the order 
of a year. Thus, the smaller the time scale r, the more precise is the normal approximation for 
rj(t, r) assuming a perfect Gaussian statistics for ryj(t, r) and the smaller are high-order cumulants 
of rj(t, r), since they quantify the departure from normality. This is the classical point of view for 
the existence of non-normal returns as studied in the financial literature. This is equivalent in our 
model to setting the parameters hi = 0.5 and Af = 0. 

Non-trivial multifractal corrections arise when Xf is not negligible and thus hi departs from 0.5. 
For the statistical properties of (both "discrete" and "continuous") returns, these multifractal 
corrections to the log-normal picture are stronger at fine scales: the finer the scale, the larger the 
values of high order cumulants. With respect to time scales, this phenomenon thus acts in the way 
opposite to the previous effect of the difference between "discrete" returns rj(t, r) and "continuous" 
returns r/j(t,r). Due to multifractality, we do not expect to observe normality for returns at small 
time scales, in contradiction with the standard geometrical Brownian model. 

To sum up, due to the difference between "discrete" returns rj(t,r) and "continuous" returns 
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rji(t,T), normality is expected at small scales while deviations from normality is predicted at large 
time scales. Due to multifractality, non-normality is expected at small scale and Gaussian statistics 
are predicted at large time scales, beyond the "integral" scale. Putting the two effects together, 
non-normality must be the rule at all time scales, with a cross-over from small to large time scales, 
from non-normality induced by multifractility to non-normality induced by the difference between 
continuous and discrete returns. Thus, we expect high order return cumulants to be non zero both 
at large scales because of the log-normality and also at small scales because of multifractality of t]. 

Let us now quantify these effects. In this purpose, we compute the cumulants of the portfolio returns 
using an expansion where the small parameter is r (measured in year), which takes multifractal 
corrections into account. Prom Eq.(^7|), the cumulants C^{t) of the portfolio can be expressed in 
terms of the cumulants C^{t) of the returns for each individual asset as weighted sums over them: 

N 

C[ir) mClir), (57) 

■t=i 

N 

(r) = < Q(r) , forn > 1 . (58) 

i=l 



Using Eqs. (pO|-|56D and from the definition r(r) = — 1, after some algebra we find, to the leading 
orders in the time expansion, 



N ^ 



Cf(r) ^ J2uj^r{f,, + -af) , (59) 
i=i ^ 

N 

C[ir) ^ Y.^lalr, (60) 

i=l 

N , s 

Clkr) ^ ^^p(^3af-18afA?ln(^)j, (61) 



N 



N 



C[{r) ^ Y.^t(l(iay-l2afxyH^)), (62) 



C[ir) ^ ^u;f (^125afr^- 240afA,Vin(^)j, (63) 

C[{t) ^ Y.^f[l296al'r' + 720afxyin\^)y (64) 

where we have set af = aj^./Ti (in year~^). We can see how the multifractal corrections arise in the 
high order cumulants. Using the typical values /ij ~ iTj ~ Aj = 0.1 and Tj = 1, one can compute, 
according to the value of r, which term from the log-normal expansion or from the multifractal 
correction dominates the cumulant. Por the first two cumulants Ci and C2, it can be shown that 
multifracal corrections are always negligible. Por the third-order cumulant C3, for time scales larger 
than a few seconds, the multifractal nature of tj is not important. However, for cumulants of order 
larger than or equal to 4, the multifractal corrections are very important for time scales as large 
as several months: for instance, using the previous values for the parameters, we find that the 
multifractal term is dominating for all times r smaller than r* ~ 6 months. 



Prom the structure of the cumulants given by Eqs (51,53,^) and their generalization to higher 



order, it is easy to show by recurrence that the multifractal correction to the 2n-th order cumulant 
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behaves, for n > 1, as 

CLir) ^ E wf-ar (-1)-^ In-^ ( J j , (65) 

to leading order. As a consequence, the corrections to the Gaussian description brought by the 
multifractal model are all the more important than the order of the cumulant is higher, i.e. than 
one looks further in the tail of the return distribution. We retrieve a common observation that the 
deviations from the Gaussian model are all the more important for large risks quantified by the 
behavior of the tails of the portfolio return distribution [Sornette, 1998; Sornette et al, 2000]. 

If we neglect the higher order cumulants in the previous expansion, we retrieve that the portfolio 
return Cf(r) and variance C^(r) are both proportional to time scale r. Hence, the portfolio 
optimization of the asset weights Wi are independent of the investment horizon r in the absence 
of the riskless asset. This is the well-known result underlying Gaussian portfolio optimization 
[Markovitz, 1959; Merton, 1990]. In contrast, in the hierarchical model, cumulants of order 4 and 
higher have a different dependence on the investment time horizon r. This implies that large risks 
exhibit a non-trivial time dependence: it will therefore not be possible to optimize the asset weights 
Wi for all time scales r simultaneously. This is the novel ingredient captured by our hierarchical 
model: a portfolio optimization which is concerned with large risks has to optimize its investment 
time horizon in a manner that we are now going to investigate. 



3.3 Excess kurtosis 

The absence of volatility correlations and the validity of the log-normal model would imply, as 



already said, that the variance of the volatilities vanishes. As seen from expressions (62-p4|), all 



cumulants of order n larger than 2 would then be of order cr^^" ^V""^, i.e. very small at small 
time scales, as expected for a Gaussian distribution of returns. Actually, for instance for the US 
S&P500 index, there is a residual hierarchical correlation structure, quantified hy h = 0.518 and 
A^ = 0.018 and T ~ 1 year. With these values, we predict an excess kurtosis by using Eq. ( |60D and 
I) for a single asset (the S&PSOO index) 



I60-V - 12A2 ln(T) . (66) 

2 



This expression makes clear that the size of the variance A^ 7^ of the volatilities at different 
scales quantifies the distance from the log-normal paradigm. The excess kurtosis at small time 
scale is negligible for A^ = 0. In contrast, in the presence of multifractality (A^ 7^ 0), the excess 
kurtosis k is predicted to be very large at small time scales. The dependence of the excess kurtosis 
is shown in Fig. |5[ This figure shows that the excess kurtosis given by Eq. ( |66|) can be represented 
approximately as a power law decay k ~ (r/T)^'^'^ over more than three decade with a small 
exponent ~ 0.2, in excellent agreement with previous determinations [Dacorogna et al, 1993; Ding 
et al, 1993; Bouchaud et al., 2000]. To confirm the relevance of Eq. (|66|) , one thus needs to 
investigate the excess kurtosis for time scales large enough as compared to the integral scale T. It 
is important to contrast this decay of the excess kurtosis k as a function of the time scale r with 
what one would expect from a model without correlations across scales: in that case, all cumulants 
are linear in r and k oc 1/r. The anomalous law (]66| ) is a clear signature of long-range correlations 
in the volatilities. 
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4 Portfolio optimization with time-scale dependent risks 



4.1 The expected utility approach 

Starting the period with initial capital Wq > 0, the investor is assumed to have preferences that are 
rational in the von-Neumann-Morgenstern [1944] sense with respect to the end-of-period distribu- 
tion of wealth Wq + 5S. His preferences are therefore representable by a utihty function u{Wo + dS) 
determined by the wealth variation 6S at the end-of-period r. The expected utility theorem states 
that the investor's problem is to maximize E[u(Wo + 6S)], where E[a;] denotes the expectation 
operator : 

/ + 00 
d6S u{Wo + 6S) P^{6S) . (67) 
-Wo 

The utility function u{W) has a positive first derivative (wealth is prefered) and a negative second 
derivative (risk aversion). Use of the utility maximization in portfolio optimization can be found 
in [Levy and Markowitz, 1979; Kroll et al, 1984]. 

Here, we consider a simple case of a constant absolute measure of risk aversion —u"/u' = a (where 
the primes denote the derivatives), for which u{W) = — exp(— aM^). With a very good approxima- 
tion for large initial wealths, we can take Ps{5S < —Wq) ~ 0. This gives 

/ + 00 
e-''^^P^{6S)d{6S) , (68) 
-oo 

which is nothing but the Fourier transform of the probability distribution function (with imaginary 
argument k): 

E^u{Wo + SS)] = -e-'*^«Pr(^ ^ (gg) 

By definition of the cumulants, this reads 

E^u{Wo + SS)] = -e-'^^o exp (^-aCf (r) + ^ iz^Cf (r)^ . (70) 

Maximizing the expected utility E'^[?x(VFo+'5'S')] thus amounts to minimizing the argument — aCf (t)+ 
J2n=2 ^~(ny. ^ni'^) second exponential in the r.h.s. of ( [70| ) with respect to the weights Wi. 

Keeping in mind that the time r is the small parameter of the problem, we express the optimal 
asset weights Wi as equal to the optimal Markowitz results denoted w° (obtained by droping all 
multifractal corrections) plus the multifractal corrections. Keeping only the cumulants of order up 
to 4 (this approximation is valid not only in the small time limit but for a risk aversion a not too 
large). The problem is thus to minimize 

N N 1 a'^ 

i=l i=l i=l 

3 N 4 N . V 

- y E ^y^'^f + ii E - 12^'^^' J)) (^1) 

with respect to the asset weights tfj's. The term a is a Lagrange multiplier ensuring the normal- 
ization X^i^i Wi = 1. This optimization amounts to finding the roots of a third order polynomial 
and one can thus obtain closed expressions for the weights Wi. 
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In order to quantify the influence of multifractal corrections and to handle simple expressions, let 
us neglect the corrections due to the difference between "continuous" and "discrete" returns which 
are unimportant at small time scales up to 6 months. In this approximation, the problem is to 
minimize 



N N ^2 N ^/i N ^ 

- a^Wi- a^WiT^Mi+ + — ^wf T af - —'^wfl2crt\^,T'^\n{—). (72) 
1=1 i=i i=i i=i * 



Then, in absence of the last terms proportional to a^, i.e. for all A?'s equal to zero, we get the 
"Markowitz" solution 

Q afijT + a 



(73) 



where a is determined from the normalization condition. In the simple case where = cr^ is the 
same for all assets, this gives 

< = ^ + ^(^J-(/^)) ' (74) 

where (/u) is the mean return averaged over all assets: assets with better than average returns have 
thus more weight in the portfolio, with a leverage controlled by the risk cr^. 



Using this solution (|73|), we get the general solution of the weights Wi up to second order in powers 



of A 



where 



Wi = 



1 + AiX 



2 



Ai = 2a^afw^iThi{— 



(75) 



(76) 



The expression of the weigths valid to first-order in the multifractal corrections for r < r* as a 
function of time horizon r is then 



Wi = Wi 



1 + 2aV>,"A;rln(— ) 



(77) 



Returning to the general solution ([75|), we see that assets with integral time scales Tj ~ er (where 
e PS 2.718 is the base of the natural logarithm) will be the most depleted compared to the Gaussian 
solution, as the factor Ai is negative with a maximum amplitude for r/Tj = 1/e. Both for small 
r/Tj and for r/Tj approaching r*, the solution is close to the Gaussian solution as it should : small 
r/Tj do not lead to large absolute risks; r/T, r* leads to the Gaussian regime r/Tj = 1 for rji. 
The worst case occurs for intermediate values of the time horizon compared to the integral time 
scale. Such stocks will be unfavored in the portfolio selection. 

This result is actually more general than this section would lead us to believe. This can be seen 
from the structure of the cumulants of order 2n given by Eq. (|65|): the dependence in r/Tj given by 

{y-^ In (t")) shows that, for large n, the cumulants are the largest when \{t /Ti)\n{T /Ti)\ 
is the largest. This is exactly the same term that controls the corrections to the Gaussian case 
quantified by the parameters A^. 



4.2 Efficient frontiers for multi-period portfolio optimization 

As recalled in the introduction, the multiperiod portfolio problem as been addressed by several 
studies. This problem is a natural application of our multi-scale description of returns. The 
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problem we investigate is to minimize a risk measure represented by a cumulant of order 2n for 
a fixed mean return. For instance, one can choose to minimize the value of the cumulant C4 at 
fixed mean return Ci (this amounts to find the set of Wi 's that minimize the pseudo- utility function 
—aCi + = —aJ2iLi WiC\ + Y^f=i wfCl), thus defining the Ci — C4 efficient frontier [Andersen 
and Sornette, 2000]. The previous cumulant perturbative expansion at different time scales r allows 
us to estimate the shape of the generalized optimal frontiers for all periods, given a fixed horizon 

For the sake of simplicity, we will suppose that all the assets are characterized by the same integral 
time values Tj = T = 1 year, which is a reasonable value as revealed by the empirical analysis of 
section 2. Denoting the investment horizon and Np the number of periods, we will consider the 
portfolio returns at scale r = Th/Np. Moreover, we will explicitely make the so-called "rebalancing 
assumption": periodic rebalancing supposes that, at the end of each period, portfolio composition 
is adjusted in order to restore the original weights. We will assume that doing so the returns 
associated to each period are statistically independent. 

Once again, in order to quantify the influence of multifractal corrections and to handle simple ex- 
pressions, we will neglect here the corrections due to the difference between "continuous" and 
"discrete" returns which are unimportant at small time scales up to 6 months (above which 
they are dominated by the log-normal corrections). One can easily show using rj ~ rji and the 
"rebalancing assumption" that the cumulants of the portfolio returns for a given multiperiod 
strategy have exactly the same expression as in Eqs. (|59|)-(|64[) where ^^r and af^r"^ become 
respectively ^iT^ and af^ThT^"'"^^ and the multifractal corrections xf^"" ^V" ln'-"'~^'*(r) become 
TfiXf^^ ^V^"~^^ ln'-"'~"'^''(r). Using this rule, it is easy to show that, if one estimates the risk using 
the variance of portfolio returns, the results are independent of the number of periods at fixed 
horizon Th- Fig. |^ represents the standard efficient frontier Cf {T^, r) as a function of C2{Th,T) for 
various r's. As expected, it is clearly apparent that the efficient frontiers are independent of time 
r, i.e. of the number of periods. This is nothing but Tobin's result [1965] that the single-period 
minimum variance set and the multi-period minimum variance set are identical. 

In contrast, the generalized efficient frontiers with a risk measure represented by higher order 
cumulants will involve multifractal corrections and both portfolio composition and efficient frontiers 
will depend on the number of periods. Let us illustrate this results for the minimization problem. 
In that case, the optimum weight of asset i can be written for r < r* as: 

where a is the Lagrange parameter associated with the fixed mean return. For a large number of 
periods (i.e. small time scale r), assets with small multifractal parameter Xf are prefered. Using 
these weights. Fig. |^ shows the efficient frontiers Cf {Th,T) as a function of C[{Th,T) for various 
periods. We have considered the simplest case of a single risky asset with parameters /i for the 
mean return, a for the standard deviation and for the multifractal parameter. Denoting /xq the 
return of the riskless asset, the parametric equation for the efficient frontier for an horizon Th and 
a number of periods T^/t is: 

1 

Cf(n.r) = ,„T.+ (-j^[^j^)"(,-,.,)r.. (79) 
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Fig. exhibits a measurable dependence on the time scale r, in contrast with the standard mean- 
variance efficient frontier Cf{Th,T) as a function of C2 {Th,T). For a fixed risk , the return 
is seen to be increasing at both ends as a function of time scale r/T, i.e., for r/T — > and for 
r/T — > T*. There is thus a worst choice for the time- horizon and for the rebalancing of the portfolio, 
approximately given by one fifth of the integral time scale T. For this time scale, the return is 
minimum for a given (large) risk and the risk is maximum for a given return. One the other hand 
one can see that the shorter the rebalancing period, the lower is the risk for a given return. As 
a matter of facts, the risk is here quantified by the fourth cumulant which goes to zero as r ^ 
(Gaussian statistics). The limiting factor for small r's is then transaction cost. Fig. 8 makes 
more precise this effect by plotting Cf {T^, r) as a function of r/T for a given value = 2 of the 
accepted risk. 



5 Distribution of returns of a portfolio of correlated assets within 
the multifractal multivariate cascade framework 



In order to discuss portfolio optimization in presence of correlated assets, let us make the approx- 
imation (valid for r small enough) that ^/^(t) ~ ^{{t). Moreover, for the sake of simplicity, let us 
assume that the r/j's are centered, i.e., /ij = 0. Using the notations of section 2.3, the characteristic 
function P{k) of the portfolio R{t) can be written as 

d^uMN(u-hiT),A{T)) / (^^rAA^(f,D(n)5]oD(n))e*'='^•"^ (81) 

-00 ^ ^ J —oo 

where we have denoted w the vector of the weights Wi. Let us consider the orthogonal matrices Q 
and O that diagonalize respectively the matrices A (and thus also the matrix A?^ by virtue of Eq. 
(p8|), showing that Q is independent of r) and Sq: 

A(r) = Q*A(t)Q, (82) 
So = 0*oSoOo , (83) 

where A and XIq are the diagonal matrices formed of the eigenvalues of A and Sq- Let us set 
r' = OY)~^{u)r. If we denote w' the vector Ow, the previous equation can be rewritten as: 

(i^nAA^(n-/I(r),A(r)) / dV nM(r>2je^'=^'''"'^'^^ (84) 

-00 ^ ^ J~oo 

where a^, are the eigenvalues of Sq and A/i is the standard monovariate Gaussian law. After 
calculating the second integral, we obtain the following expression for the characteristic function: 

d^uAA^ (n-/i(r),A(r)) e"~^-i' ""-^-t, . (35) 



From this equation, one can immediately see that, if all the cascades are uncorrelated, i.e, A 
is diagonal, the characteristic function is the same as in the case of independent assets where 
the weights Wi have been replaced by the weights w[. This is the same situation as in classical 
(Gaussian) portfolio theory where the case of correlated assets is reduced to the uncorrelated one 
by such a simple change of variable on the weights. 

When the covariance matrix A is non trivial, the situation is more complicated. Let us consider the 
case where each matrix A^ in the block decomposition ( [2^ ) is singular, i.e, [A^Jjj = A^. This case 



18 



corresponds to the existence of L degrees of freedom in the problem for which one can compute the 
characteristic function under the form: 



1=1 



This case corresponds to L independent assets whose integral scale variances are Nia^^. 
Expression ( |86|) is the product of L terms of the form (^) already encountered for the case of a 



single asset in section 3. The moments and cumulants can thus be calculated analytically. 



6 Conclusion 



We have extended a statistical model of price returns and volatility correlations based on the idea 
of information cascades from large time scales to smaller time scales. Empirical tests performed on 
intra-day as well as daily data on the CAC40 and S&PSOO indices, on their constituting stocks as 
well as on bonds and currencies validate satisfactorily the model. The calibration give a robust and 
seemingly consistent value for the two key parameters: the integral time-scale is found in the range 
of one to two years and the variance of the multiplicative kernel is approximately 0.02 for all stocks 
and indices that have been investigated. Our results show that the evidence for multifractality is 
fully consistent with a simple cascade origin, flowing from large time-scales to shorter time-scales. 
The multifractal cascade model offers an intuitive explanation for the observation that the detection 
of "abnormal" states or crises in the stock market requires an index constructed over many different 
horizon times [Zumbach et al., 2000]. 

We have also offered an extension of the cascade model into a multi-variate framework to account for 
correlations between assets, in addition to the correlations in time-scales. Future works will exploit 
this novel formalism, in particular for multi-period portfolio characterization and optimization. 

In a second part, we have shown how to characterize the distribution of returns of a portfolio 
constituted of assets with returns distributed according to such multifractal cascade distributions. 
In particular, explicit analytical expressions for the first six cumulants are offered. We also show 
that, within a utility approach with a constant absolute measure of risk aversion (exponential util- 
ity function), the problem of portfolio optimization amounts to maximize a sum over cumulants 
weighted by powers of the risk aversion coefficient. Working in the space of (return, fourth-order 
cumulant) or of (return, sixth-order cumulant) generalizes the mean-variance approach and under- 
lines the impact of the investor horizon-time. The most important consequence of the theory is that 
the optimal portfolio depends on the time-scale. In addition, it is not possible to simultaneously 
optimise all the components of risks with respect to the choice of the investment time scale. This 
result extends to the investment horizon dimension previous results obtained from a decomposition 
of the risk into a spectrum from small to large risks quantified by the cumulants of the distribution 
of portfolio returns [Sornette et al., 2000]. 

In principle, our portfolio theory allows one to quantify how much diversification can be obtained 
by buying different assets and managing them optimally with respect to their possibly differ- 
ent integral time scales: indeed, our results suggest that reallocation of assets in the portfolio 
should be performed with different time-horizon depending upon the assets. This is related to the 



19 



"time-diversification" concept introduced by Martellini [2000] for option hedging in the presence of 
transaction costs. A quantification and tests of these strategies will be reported elsewhere. 
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Figure 1: (a) the Standard & Poor US index (S&;P500) time series (left panel) and the Japanese 
Yen time series (right panel) with a sampling time of ten minutes are shown as a function of 
time over more than one (resp. two) decades. Intraday "seasonal" effects (large volatility at the 
open and close) are removed before calculating the moments, (b) log2 (M(g,T)) versus log2(r) for 
q = 1,2,3,4,5. The slope of these curves provides an estimate of Cio.)- (c) log2(Af(g, T)/Af''(l, r)) 
versus log2(r) for q = 2,3,4,5 (top to bottom); deviation from a constant is the signature of 
"multifractality" . (d) Cil) spectra for the two analyzed series. As illustrated in the insets where the 
linear part of the spectra has been removed, they arc almost perfectly fitted by a pure parabola with 
quadratic coefficient = 0.018 (S&P500) and = 0.032 (Japanese Yen). We find h-\^ = 1/2 
to a very good approximation (see text for an explanation). 
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Figure 2: This figure shows the autocorrelation function defined by Eq. (^Oj) of the "magnitudes" 
LOi{t,T) = ln(|r7j(t, r)|) for r = 10 minutes for the SP&500 (o) and Japanese Yen (x). A hnear 
regression of these functions versus the logarithm of time lag IS.t enables us to evaluate A? and Tj 
for each time series. 
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Figure 3: Panel (a) represents in dotted line the magnitude correlation function for the pair of 
stocks CCF/MICHELIN of the French CAC40 index. The mean of such correlation functions over 
all stock pairs of the CAC40 index is the continuous line. The slope and the intercept of these 
functions allow us to evaluate \fj and Tij. Panel (b) shows the histogram of Xfj for all index pairs. 
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Figure 4: ln{Tij) versus Xfj : each circle represents a pair of stocks taken from the French 
CAC40 index (resp. from the American Dow Jones index for the darker circles). For both indexes, 
the correlation time-scales and multifractal coefficients are approximately clustered on the same 
~ 0.02 and T ~ 1 - 2 years. 
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ln(T/T) 

Figure 5: Excess kurtosis k given by the prediction (^) of the multifractal cascade model with 
parameters = 0.018 and cr^ = 1.8 10""^ (see S&PSOO data in Table 1) in logarithmic scale as 
a function of t/T also in logarithmic scale. The decay is not a pure straight line but can be 
approximated by one with an average exponent close to 0.2 for small r's as previously observed 
[Dacorogna et al, 1993; Ding et al, 1993; Bouchaud et al., 2000]. 
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Figure 6: Standard efficient frontier for multiperiod portfolio: Ci{Th, r) as a function of C2{T}^, t) 
for various time-scales r's. Units are arbitrary. The shape of the frontier illustrates the indepen- 
dence with respect to the time horizon. 
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Figure 7: Generalized multiperiod efficient frontiers Ci{Th,T) as a function of C[{Th,T) for 
various r/T's in the case of a single risky asset. The return of the unrisky asset has be chosen to be 
/io = 0.05 and the parameters of the risky asset are fi = 0.15, a = 0.1, T = 1 year and = 0.02. 
T/i is set to 1 year. 
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Figure 8: Return quantified by C[ {T^, r) as a function of the reduced time scale t/T for a given 
value of C[ = 2. The data are the same as in Fig. 7. One can note that for a given risk, return 
is better at very small r's than at very large ones. Moreover, there is a worst horizon for which 
return is minimum. 
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